Energy levels and far-infrared spectroscopy for two electrons in a 

semiconductor nanoring 



Hui Hu* 

Department of Physics, Tsinghua University, Beijing 100084, P- R- China 

Jia-Lin Zhu 

Department of Physics, Tsinghua University, Beijing 100084, P- R- China 
Center for Advanced Study, Tsinghua University, Beijing 100084, P- R- China 

Jia-Jiong Xiong 

Department of Physics, Tsinghua University, Beijing 100084, P- R- China 

Abstract 

The effects of electron-electron interaction of a two-electron nanoring on the 
energy levels and far-infrared spectroscopy have been investigated based on a 
model calculation which is performed within the exactly numerical diagonal- 
ization. It is found that the interaction changes the energy spectra dramat- 
ically, and also shows significant influence on the far-infrared spectroscopy. 
The crossings between the lowest spin-singlet and triplet states induced by 
the coulomb interaction are clearly revealed. Our results are related to the 
experiment recently carried out by A. Lorke et al. [Phys. Rev. Lett. 84, 2223 
(2000)]. 

PACS number(s): 71.15.-m, 73.61.-r, 78.30.-j, 73.20.Dx 



* Email: huhui@Phys .Tsinghua.edu . cn 







I. INTRODUCTION 



Rapid progress in nanostructure technology has made it possible to fabricate various 
nanometer quantum devices which have potential applications in microelectronics. Such 
ultrasmall devices contain only a few electrons and the electron-electron interaction is pro- 
posed to be of great importance to theirs energy level structures and optical properties. It 
leads to a number of new quantum phenomena. One of the most interesting phenomena is 
the spin oscillation of the ground state in an external magnetic field, which is due to the 
interplay between three energies: the confinement potential, the Zeeman energy, and the 
electron-electron interaction. The simplest case, a pillar few electrons quantum dot (QD) 
with a parabolic potential, has been extensively investigated, and the spin oscillation in a 
QD is well understood.S^i The theoretical predictions are also confirmed in the experiment 
through the conductance measurements in the finite drain-source voltage regime.@ 

The semiconductor quantum ring is another interesting example. In 1993, D. Mailly et al. 
measured the persistent currents in a mesoscopic single GaAs ring induced by a magnetic flux 
threading the interior of ringjl and the experimental results have attracted a lot of theoretical 
interests. One of the basic questions addressed by many theoretical explanation is concerned 
with the role of the electron-electron interactioniHil For a narrow-width rings, an adiabatic 
approximation allows one to decouple the radial motion from the angular motions and to 
arrive at analytical solutions for the wave functions and the energy spectra.00 As shown 
by L. Wendler et al. ,00 the interplay of the Coulomb repulsion between the electrons and 
the confining potential forms a relatively rigid rotator with internal azimuthal excitations 
and confined radial motions, i.e., the picture of a rotating Wigner molecule. However, at the 
moment, the influence of the electron-electron interaction in a finite-width and nanoscopic 
quantum ring, i.e., nanoring, is still less well understood.llHll 

Very recently, using the self-assembly techniques, A. Lorke and collaborators demonstrate 
the realization of nanoscopic semiconductor quantum rings inside a completed field-effect 
transistor (FET) structure.EH^I Quite different from the conventional sub-micron meso- 
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scopic structures, the nanorings are in the true quantum limit. By applying two comple- 
mentary spectroscopic techniques, capacitance-voltage (CV) spectra and far-infrared (FIR) 
spectroscopy, they investigate both the ground state transition and excitation's properties 
of these two-electron nanorings in a magnetic field perpendicular to the plane of rings. Al- 
though the main experimental results can be qualitatively explained by the single electron 
picture, some contradictions remain, i.e. the coulomb interaction energy estimated roughly 
20 meV, is too large to be ignored safely. Hence, more in-depth theoretical works are de- 
sirable, especially in view of the existences of the very strong electron-electron coulomb 
interaction. 

In this paper, we would like to study the energy levels and FIR spectroscopy of a two- 
electron nanoring, and pay special attention to the effects of the Coulomb interaction. First 
of all, for a ring-like confinement potential, the total Hamiltonian cannot be separated into 
the center-of-mass and the relative-motion terms. We develope a new theoretical method to 
handle this non-separability. It consists of the well-known series solution method,§'0 which 
is effective to solve the single particle problem, and the exact diagonalization method.il'0 On 
the other hand, we show that the electron-electron interaction can change the energy levels 
significantly. An obvious feature induced by the interaction is the intersection between the 
lower levels. It presents the spin oscillation of the ground state in an external field. Moreover, 
the results obtained by A. Lorke et al. can be explained more realistically in our model. 
Further, the generalized Kohn theorem breaks down in the FIR spectroscopy due to the 
mixing between the center-of-mass and the relative-motion modes as mentioned above.@il 
We outline here two prominent features of the FIR spectroscopy caused by electron-electron 
interaction: splitting and dis-continuous drops of the resonance energies, which are suggested 
to be detectable with circularly polarized far-infrared light. 

The organization of this paper is as follows. In the next section, we introduce the model 
and Hamiltonian briefly. In Sec III, the main procedures of the numerical calculation are 
outlined. Sec IV is devoted to the results and discussions, followed by a summary in Sec V. 
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II. THE MODEL AND HAMILTONIAN 



The nanoring is considered to have two electrons with an effective conduction-band-edge 
mass m* moving in the x — y plane, and a ring-like confining potential can be introduced 
U(r) = ^m* e uj\ [r — R f^E where ujq is the characteristic frequency of the radial con- 
finement and Ro is the ring's radius. This system is subjected to a perpendicular uniform 
magnetic field which is described by a vector potential A (r) = \B x r in the symmetric (or 
circular) gauge. The resulting Hamiltonian is given by 

where = (xj, y^) and Pi = —ih\7i are respectively the position vector and momentum 
operator of the 2-th electron with charge — e. Eq is the vacuum permittivity and e r is the 
static dielectric constant of the host semiconductor. In addition, there also exist the spin 
interaction term with the magnetic field 7i sp in = g* \lb \ S\ + S2 ) ' ^1 wnere 9* is the Lande 
factor, and /j,b is the Bohr magneton. Due to the small g* in semiconductor, in general, 
H, S pi n is very small and can be ingored safely (i.e., for g* = —0.44 in GaAs materials and 
B — 10 T, the typical value of Tispin is 0.25 meV). It is worthy to note that in the limit 
case Rq = 0, the nanoring simply reduces to a parabolic quantum dot, which can be solved 
trivially by separating the Hamiltonian into center-of-mass and relative-motion terms. The 
appearance of Ro breaks this separability, and it becomes more complex. 

We apply the exact diagonalization method by constructing the basis with the single 
particle wavefunctions of Hamiltonian 

H s = ^(p + eA(f)y + U(r). (2) 

These wavefunctions labelled by the radial quantum number n and orbital angular- 
momentum quantum number m have the form 

ipnm(r) = R nm (r) exp (imip) n = 0,1,2..., m = 0,±l,..., (3) 

where the radial part R nm will be solved exactly by using the series expansion method.! 
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III. FORMULA AND CALCULATION METHODS 



A. Series solution 



For the sake of convenience, we use the effective atomic units, in which the effective 
Rydberg R* = nb0 ™ ee — and the effective Bohr radius a* = 47r£ ° £ ' 2 fi2 are taken to be the 
energy and length units, respectively. Then, the Hamiltonian (|J) has the form 

H, = -V 2 + ~ 7fe V + -^l (r - R ) 2 + lb L z , (4) 

where the magnetic field % is measured in the unit with cyclotron frequency oj, - ' h 



*\l/2 



7;,Z/ 2 is the Zeeman term, and r y d 1//2 = (^j^j is related to confinement region of the 
electrons. It is interesting to note how large the units of semiconductor materials are. For 
GaAs materials, for example, R* = 5.8 meV, a* = 10 nm, and 7& = 1 corresponds to 
B = 6.75 T. 

Now, we have to solve the Schrodinger-like equation 

T-L s [Rnm{r) exp (irrvp)] = E nm [R nm {r) exp (irrup)] (5) 

to obtain the energy E nm and radial part of wavefunction R nm (r). It is easy to find the 
equation satisfied by the function R nm (r) 



d 2 1 d 
dr 2 r dr 



Rnm (r) = 0. (6) 



We are prevented from analytically exact solutions of the eigenvalue problem because Eq. 
(g) with suitable boundary conditions is beyond the analytical problem of confluent hyper- 
geometric equations. However, we can use the method of series expansion to obtain the 
series forms in different regions of Eq. (|J) and the exact values of i? nm .i'0 

It should be noted that r = and r = +oo are respectively the regular and irregular 
singularity points of Eq. @. So it is natural to divide the whole region [0, +oo) into three 
parts, [0, +oo) = [0, r ) U [r , Too) U [roo, +oo), where r and r ro are the two dividing points. 
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With suitable adjustment of r$ and r^, very high numerical precision can be archived. In 
all the three kinds of regions, the function R nm (r) is found to be the form: 



Rnm (j) * 



oo 

Ar l «n r ™ ^ r < r , 

n=0 

oo oo 

Ct E Ci„ (r - + A £ d m (r - Ri^r< R i+1} (7) 

n=0 n=0 

oo 

r^a € r < +00, 



/ 2 \ 00 

B exp ( — ^r 2 + t^RqT 1 r s £ & n ^~ n 

^ " ' n=0 



where 7 = V7b~+7l; I = \ m \ an d s = <yEnm ~ mib ^ — 1 — ^t-Rq. In order to improve precision, 
we further divide the region [r , r^) into N pieces, denoted by R4 (i — 1, N), here i?i = r 
and i?Ar = Too. A, Ci, Di and .B are constants, a n , Ci n , di n and h n are expanding coefficients 
and can be determined by the recurrence relation coming from Eq. (||), the initial values of 
coefficients are chosen to be do = 1, Qo = da = 1, en = d i0 = and 60 = 1- 

Using the matching conditions at r = R$ (i — 1, iV), and the 2x2 transfer matrices, 
we can deduce the equation for eigenenergies E nm easily. Then the constant A, Ci, Di and 
B can be evaluated by normalization condition, and Rnm (r) is obtained numerically. 

To close this subsection, it is interesting to point out that the method mentioned above 
is very suitable for numerical calculations, and can be modified to handle various differential 
equations similar to Eq. (0). 



B. Exact diagonalization 

Once the single particle wavefunctions are obtained, we go ahead to construct the basis 
of the two-electron wavefunctions. It is obvious that the total z component of the angular 
momentum operator of two electrons, L Zyto tai = L z \ + L z2 , is a constant of motion, i.e., 
Lz,totai — is valid, from which the rotational invariance of the problem follows. Thus, 
one can work in one subspace labelled by a good quantum number L instead of the whole huge 
Hilbert space. We denote the corresponding two-electron wavefunction by ipL (fici; r^a-i) 
which depends on the spatial coordinates {r*j} and the spin coordinates {<jj} (i = 1,2). 
Because the Hamiltonian of Eq. (P does not depend on the spin operator, in our two 



electrons case, the wavefunction can separate into the orbital part ipL (ri,f 2 ) and the spin 
part xs (0i,<r 2 ): 



fpLS (^>ii r 2 a 2 ) = ipl lX2 (fi, r 2 ) Xs (d, ^2) 



(8) 



with 



^l iA2 (n, r 2 ) = c Vmmi (n) Vw 2 (r 2 ) + (-)° ^„ 2m2 (ri) ^ mmi (f 2 ) 



(9) 



where the normalized wavefunction ipis r 2 a 2 ) is further labelled by the total spin 
number S — or 1, corresponding to the singlet and triplet state. Aj (z = 1,2) stands for 
the quantum number pair (n^mj) . L = m\ +m 2 . ip nm (r) is the single particle wavefunction. 

1 as (riimi) 7^ (n 2 m 2 ) , and c = | as (r^mi) = (n 2 m 2 ) . Obviously, the wavefunction 
constructed above satisfies the antisymmetric condition, 



VuiPls (rioi; f 2 cx 2 ) = -^ L5 (fto; r 2 a 2 ) , 



(10) 



where 7 3 i 2 is the permutation operator. 

In the next step, we diagonalize the Hamiltonian ([]]) numerically in a restricted config- 
uration space. The space is constructed by choosing the wavefunctions having the form (||) 
in the lowest / levels. The secular equation of finite degree / is given by 



det 



eP - E ) S i:j + Aij 



0. 



i,j = 1, ■-.,/, 



(11) 



where Ef' = E n i m i + E n i m i is the single particle energy, i = (n\m\; n l 2 m l 2 ) and j = 
(n{m{] n^rn-^ represent the quantum number levels. is the matrix element of electron- 
electron interaction in unit of R*. 



\r 1 - r 2 \ 



-f A 



(12) 



where 
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dfidf 2 ^* n i mi {ri)V n i m i (r 2 ) 



(13) 



and 

A tj = J J ^l^Cimf fa) V£*m* ^ |^~^[ ^^ fa) ' ( 14 ) 

A?-, A?- can be computed numerically from the series solution of ip nm (r). 

By diagonalizing the secular equation ([11]) in each subspace (S,L), we obtain the m-th 
energy level E m and the corresponding two-electron wave-function: 

$ m (fian r> 2 ) = ^ ^ lA2 ^ lAa (fx, F 2 ) X 5 K «r 2 ) . (15) 

Ai A2 



C. Optical absorption 



In the electronic dipole approximation,Q the absorption coefficient is given by 



aM = cu^2 \\f 

fi 



e ■ d 



(16) 



where e is the complex polarization vector of the spatially constant external electronic field, 
d = — (fi + F2) is the electronic dipole operator of two electrons, and is the equilib- 
rium Fermi distribution. The summation is over all the two-electron eigenstates, LOfi is 
proportional to the energy difference between the initial state | i) and final state | /) : 
ujfi = u>f — uji = (j-^fi— - J ■ c is a constant factor. Restricting ourselves to zero temperature 
T = OK, Eq. (TT§) reduces to 
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alu) = cjjj 



e ■ d 







S(u- Ufa) , 



(17) 



where | 0) and | /) represent the ground state and excited state respectively. 

For circularly polarized light we have e = « (1, ±i) , from which it follows e ■ d 



-\ \{r x e ±l ^ + (1 -> 2)). Using Eq. (pi), we obtain 



+ (1-2) , (18) 



A?A^ a{a| 



where the single particle matrix elements g^ a , are defined by 



^AA' = l r ex P (± i( f)\ = $m±l,m> I r 2 Ry (r) i? A M t&" 



(19) 



and | A) represents the single particle wave-function ip nm (r). It is obvious that the absorption 
coefficient satisfies the dipole selection rule AL = ±1. Substituting Eq. (ITS) into Eq. ( |TTD 
and taking 



8(u- Uf ) 



T/n 



(uj - Ufo) + r 2 

where T is a phenomenological broadening parameter, we arrive at 



(20) 



UJ 



E 



A i A 2 A{ A^ 



2) 



r/vr 



^ /0 )' + r 2 



(21) 



where ± corresponds to right and left circularly polarized lights, respectively. 

To check the numerical accuracy of the calculations we have used the /-sum rules for 
the dipole operators, which can be expressd in terms of ground-state quantitiesSil 



(a + {u) + a~ {u))duj= (0| 



e • d 



TiAe-d 



|0) = JV 



(22) 



in the effective atomic units, where c is taken to be 1, and iV = 2 is the number of electrons. 

It is important to point out that in our case, due to non-separability of the center-of- 
mass and the relative-motion modes, the generalized Kohn theorem, which means that FIR 
can only be used to excite the center-of-mass modes of electrons parabolically confined in 
circular quantum dot, will not be held further. Thus we expect our FIR absorption result 
may reflect an excitation of the relative motion of two electrons. 



IV. RESULTS AND DISCUSSIONS 



To explain the experimental measurements,^ we have taken the material parameters 
e r = 12.4 and m* = 0.067m e for GaAs in our calculations. The radial confinement strength 
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hujQ and the ring's radius Rq are chosen to be 12 meV and 14 nm respectively. The cor- 
responding width of ring is about 15 nm, which means that the electrons are confined in 
a wide ring. Thus in contrast to the rotating Wigner molecule picture in a narrow-width 
quantum ring,00 the more pronounced energy spectra and optical properites are expected. 
For the calculation in each (S, L) subspace, we first solve the single particle problem and 
save several hundreds single particle states, then pick up the suitable single particle states 
to construct thousands of two-electron states, among which only the lowest / energy levels 
are selected. Here we simply point out that our numerical diagonalization scheme is very 
efficient and essentially exact in the sense that the accuracy can be improved as desired by 
increasing /. For instance, for the ground state in (0, 0) subspace, the use of 64 basis states 
allows the precision to be within the relative convergence of ~ 10~ 4 . On the other hand, by 
checking the /-sum rule, we find that the relative error of our FIR calculation is less than 
1%. 

A. Spin oscillation 

The energy levels of two electrons in a nanoring as a function of the magnetic field have 
been plotted in Fig.l. As mentioned above, only the total spin and total angular momentum 
are conserved in our model, and then for the sake of clearness, we only plot the lowest level 
in each (S, L) subspace. In order to understand the role of the electron-electron interaction 
in the two-electron spectra better, we first describe the characteristics of the energy levels of 
two electrons in the nanoring without interaction. As shown in Fig. 1(a), for small magnetic 
field, the spectra have the characteristics of disk-like quantum dot, since the electrons are 
confined in a wide nanoring whose radius is comparable to its effective width. As the 
magnetic field increases from zero, there are minima for the states with negative angular 
momentum L. These minima are caused by the interplay between the Zeeman term and the 
ring-like confinement potential. Moreover, the level with L = increases monotonically, and 
changes more dramatically than the others. This leads to a ground state transition from 
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L = to L = —2 around B = 8 T, which also reflects the fact that the nanoring becomes 
more and more narrow with increasing the magnetic field. In the sufficient high field regime, 
it is quite safe to speculate that the levels will be the same as the levels in one-dimensional 
ring under a uniform magnetic field. 

The electron-electron interaction can significantly change the characteristics of the spec- 
tra described above. As shown in Fig. 1(b), clearly visible is the ground state transitions 
around B = 3, 7 and 10 T. These transitions are quite different from that mentioned in 
the previous paragraph: not only the total angular momentum but also the total spin are 
changed. They present the spin-singlet-spin-triplet oscillation of the ground state in the mag- 
netic field, i.e., (0, 0) — > (1, — 1) — > (0, —2) — > (1, —3) states and so on. This phenomenon is 
indeed qualitatively similar to that seen in a quantum dot. 11 

For a better understanding of the singlet-triplet oscillation, it is interesting to study the 
electron-electron interaction energies E r , defined by the difference between the energy E with 
interaction and E^ without interaction. In Fig. 2, the E r for different states are plotted as 
a function of the magnetic field. It is readily seen that the E r increase with the magnetic 
field and the ordering is as follows: E r (0,—1) > E r (0,0) > E r (0,~3) > E r (l,—2) > 
E r (l,0) > E r (l,—1) > E r (0,-2) > E r (1, — 3) .... Compared with the corresponding 
results in quantum dot,! we find that the ordering depends significantly on the form of the 
confinement potential, e.g., the ordering E r (0,0) > E r (0, —1) in quantum dot is reversed 
in our case. However, the ordering E r (0, 0) > E r (1, —1) > E r (0, —2) > E r (1, —3) is still 
preserved, thus the trivial crossover around B = 8 T in Fig. 1(a) moves to the position 
B = 3 T in Fig. 1(b). 

We comment on some other effects caused by the interaction, (z) When the electron- 
electron interaction is excluded, the states (0, L) and (1,L) (L = —1, —3) are degenerate in 
the whole regime of the magnetic field. This degeneracy can be understood by the fact that 
those states are constructed with the same single particle states, the only difference between 
them is spin part which have no influence to energy levels in the absence of interaction. 
When the electron-electron interaction turns on, the whole energy levels are shifted to high 
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energies due to the repulsive coulomb interaction and the degenerated levels split, (ii) In the 
high magnetic field regime, the energy levels appear to be more separated in the presence 
of interaction, and it indicates that the interaction effect becomes larger in this regime. 

B. FIR spectroscopy 

In the previous subsection, we have demonstrated that the electron-electron interaction 
can induce some transitions of the ground state. It may be possible to observe these tran- 
sitions in the lower energy optical absorption. Now, we present the FIR spectroscopy to 
elucidate such transitions. Here, the phenomenological broadening parameter T is assumed 
to be 0.5 meV. 

In Fig. 3, we plot the FIR absorption for the circularly polarized lights as a function of the 
magnetic field. Compared with the no interaction case, where the right and left circularly po- 
larized (labeled by a±) resonance energies are roughly given by u± = \ {\f^ + 4cJq ± cj c jin 
the low field (as shown in Fig. 4, a + and <r_ are represented by the square and filled circle 
symbols, respectively), we find some unusual features in the presence of interaction, (i) As 
the magnetic field increases, the lowest right circularly polarized <r + resonance energy shows 
discontinuous drops around 3, 7 and 10 T (see Fig. 3a). In contrast, this behavior is only 
found at B « 8 T in the case without the Coulomb interaction (see Fig. 4). Naturally, the 
discontinuous drops are originated from the spin oscillation of the ground state as mentioned 
above. Note that they are less reflected in the <r_ case, (ii) The absorption peak around 20 
meV splits into three or more subpeaks, especially in the er_ polarization case. The striking 
behavior can be understood from the mixing of the center-of-mass and the relative-motion 
modes. Certainly, it has to keep in mind that due to the mixing there are many irregularly 
spaced and near-degenerated energy levels hybridized from the relative-motion mode in each 
(S, L) subspace. The splitting can be identified as a transition from the ground state to those 
"hybrid" states, showing its two-electron characteristic. (Hi) For both circular polarizations, 
the lowest resonance energy shows slight blue shift due to the Coulomb interaction. 
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Motivated by the above prominent features, we suggest that the spin transitions of 
the ground state and the other effects induced by the electron-electron interaction can be 
observed with circularly polarized light. However, as pointed out by A. Lorke, in theirs 
experiment, the single particle states are quite accurate basis for the description of the 
measured FIR resonance. It seems in contrast to ours expectation. We argue here that 
the controversy comes from the resolution of the experiment. With low resolution, one is 
difficult to distinguish the details of the FIR spectroscopy discussed above, and only the 
profile is explored. As shown in Figs. (3a) and (3b), we can indeed observe that the profile 
of the FIR resonance can be described by the single particle picture. Therefore, we expect 
that our predictions about the electron-electron interaction may be confirmed by a high 
resolution experiment in the future. 



V. SUMMARY 

In conclusion, we have investigated the energy levels and far-infrared spectroscopy of a 
two-electron nanoring in a magnetic field. Because of the electron-electron interaction as 
well as the interplay between the magnetic field and ring-like confinement potential, the 
nanoring exhibits rich electronic structures. An obvious feature induced by the interaction 
is the intersection between the lower levels. It presents the spin oscillation of the ground 
state on the magnetic field, i.e., (0,0) — > (1,-1) — ► (0,-2) — > (1,-3) states and so on. 
This phenomenon is indeed qualitatively similar to that seen in a quantum dot, suggesting 
its intrinsic nature of zero-dimensional quantum structures. 

The profile of the FIR spectroscopy is roughly captured by the single particle picture 
as indicated by a recent experiment .0 However, the ring-like confinement potential doesn't 
allow the application of the generalized Kohn theorem. Thus the single particle picture is 
inadequate for seeing any effect due to electron-electron interaction. We have outlined here 
two prominent features of the FIR spectroscopy caused by electron-electron interaction: 
splitting and dis-continuous drops of the resonance energies. We suggest that those two 
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features can be detectable by using the circularly polarized far-infrared light with high 
experimental resolution. 
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Figures Captions 



Fig. 1. The energy levels of two electrons in a nanoring are plotted as a function of the 
magnetic field with Hujq = 12 meV and Rq = 14 nm in the (a) absence and (b) presence 
of electron-electron interaction. The spin singlet and triplet states are labelled by solid 
and dashed lines, respectively. The quantum number of each state (total spin, total angular 
momentum) are also indicated. Note that only the lowest energy level of each (S, L) subspace 
is selected. 

Fig. 2. The net Coulomb energies are plotted as a function of the magnetic field. The 
other parameters are the same as in Fig. 1. 

Fig. 3. A logarithmic 3D plot of the far-infrared absorption coefficient as a function of 
the magnetic field for (a) right and (b) left circularly polarized light. 

Fig. 4. Far-infrared absorption resonance energies in the absence of electron-electron 
interaction as a function of the magnetic field. The square and filled circle symbols are 
corresponding to right and left circularly polarized lights, respectively. 
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